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We calculate the Casimir force and free energy for plane 
metallic mirrors at non-zero temperature. Numerical eval- 
uations are given with temperature and conductivity effects 
treated simultaneously. The results are compared with the 
approximation where both effects are treated independently 
and the corrections simply multiplied. The deviation between 
the exact and approximated results takes the form of a tem- 
perature dependent function for which an analytical expres- 
sion is given. The knowledge of this function allows simple 
and accurate estimations at the % level. 



I. INTRODUCTION 

The Casimir force |lj has been observed in a num- 
ber of 'historic' experiments It has been mea- 
sured recently with an improved experimental precision 
. This should allow for an accurate comparison with 
the predictions of Quantum Field Theory, provided that 
these predictions account for the differences between real 
experiments and the idealized Casimir situation. In par- 
ticular, experiments are performed at room temperature 
between metallic mirrors and not at zero temperature 
between perfect reflectors. The theoretical expectations 
should be computed with a high accuracy if the aim is to 
test agreement between theory and experiment at, say, 
the 1% level. The efforts for accuracy are also worth for 
making it possible to control the effect of Casimir force 
when studying small short range forces |9j-|TH]. 

The influence of thermal field fluctuations on the 
Casimir force are known to become important for dis- 
tances of the order of a typical length [|F2Ht5[ 
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When evaluated at room temperature, this length At is 
approximately 7/xm. In contrast, the finite conductivity 
of metals has an appreciable effect for distances smaller 
than or of the order of the plasma wavelength Ap deter- 
mined by the plasma frequency w P of the metal (see Q 
and references therein) 
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For metals used in the recent experiments, this wave- 
length lies in the range 0.1/im-0.2/im. This means that 
conductivity and thermal corrections to the Casimir 
force are important in quite different distance ranges. 
Thermal corrections are usually ignored in the sub-/xm 
range where the effect of imperfect reflection is signifi- 
cant whereas the conductivity correction is unimportant 
above a few fim where the effect of temperature becomes 
appreciable. This explains why the 2 corrections are usu- 
ally treated independently from each other. When an 
accurate comparison between experimental and theoret- 
ical values of the Casimir force is aimed at, the error 
induced by this approximation has however to be pre- 
cisely evaluated. Furthermore, the region of overlap of 
the two corrections is precisely in the /xm range, which 
is a crucial distance range for the comparison between 
experiment and theory. 

The purpose of this paper is to give an accurate evalua- 
tion of the Casimir force F taking into account finite con- 
ductivity and temperature corrections at the same time. 
To characterize the whole correction, we will compute the 
factor rjF describing the combined effect of conductivity 
and temperature 
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i^Cas is the ideal Casimir force corresponding to perfect 
mirrors in vacuum. L is the distance between the mir- 
rors, A their surface and h and c respectively the Planck 
constant and the speed of light. We will also evaluate the 
factors associated with each effect taken separately from 
each other 
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F p is the Casimir force evaluated by accounting for finite 
conductivity of the metals but assuming zero tempera- 
ture and F T is the Casimir force evaluated at temper- 
ature T for perfect reflectors. Of course r]^ depends on 



the ratio and rj^ on the ratio 



Now the question raised in the previous paragraphs 
may be stated precisely: to which level of accuracy can 
the complete correction factor rfp be approximated as 
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the product of the factors rjp and 77J ? To answer this 
question we will evaluate the quantity 
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which measures the degree of validity of the approxima- 
tion where both effects are evaluated independently from 
each other. We will give an analytical estimation of this 
deviation which may thus be taken into account with- 
out any difficulty. We will also give the same results for 
the Casimir energy by defining a factor 77^ measuring the 
whole correction of Casimir energy due to conductivity 
and temperature and then discussing the factors r^, and 
and the deviation Sb in the same manner as for the 
force. 

Some additional remarks have to be made at this point. 
First, recent experiments are not performed in the plane- 
plane but in the plane-sphere configuration. The Casimir 
force in this geometry is usually estimated from the prox- 
imity theorem p7|-|2l|. Basically this amounts to evalu- 
ating the force by adding the contributions of various dis- 
tances as if they were independent. In the plane-sphere 
geometry the force evaluated in this manner turns out to 
be given by the Casimir energy evaluated in the plane- 
plane configuration for the distance L being defined as 
the distance of closest approach in the plane-sphere ge- 
ometry. Hence, the factor ?7e evaluated in this paper for 
energy can be used to infer the factor for the force mea- 
sured in the plane-sphere geometry. Then, surface rough- 
ness corrections will not be considered in the present pa- 
per. Finally the dielectric response of the metallic mir- 
rors will be described by a plasma model. This model 
is known to describe correctly the Casimir force in the 
long distance range which is relevant for the study of 
temperature effects. Keeping these remarks in mind, our 
results will provide one with an accurate evaluation of 
the Casimir force in the whole range of experimentally 
explored distances. 



II. CASIMIR FORCE AND FREE ENERGY 

When real mirrors are characterized by frequency de- 
pendent reflection coefficients, the Casimir force is ob- 
tained as an integral over frequencies and wavevectors 
associated with vacuum and thermal fluctuations p2[ . 
The Casimir force is a sum of two parts corresponding 
to the 2 field polarizations with the two parts having the 
same form in terms of the corresponding reflection coef- 
ficients 
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rj_ (respectively rn) denotes the amplitude reflection co- 
efficient for the orthogonal (respectively parallel) polar- 
ization of one of the two mirrors. The mirrors are here 
supposed to be identical, otherwise r\_ should be replaced 
by the product of the two coefficients, u is the frequency 
and k the wavevector along the longitudinal direction of 
the cavity formed by the 2 mirrors. F [u] is defined for 
positive frequencies and extended to negative ones by 
parity. 

The Casimir force (0) may also be rewritten after a 
Fourier transformation, as a consequence of Poisson for- 
mula in 
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The contribution of vacuum fluctuations, that is also the 
limit of a null temperature {ut — > 0) in (^), corresponds 
to the contribution m — in (R) 
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Hence, the whole force (0) is the sum of this vacuum 
contribution m = and of thermal contributions m =/= 0. 

We will consider metallic mirrors with the dielectric 
function e {iu) for imaginary frequencies given by the 
plasma model 
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cjp is the plasma frequency related to the plasma wave- 
length Ap by (J2j> . For the metals used in recent experi- 
ments, the values chosen for the plasma wavelength Ap 
will be 107nm for Al and I36nm for Cu and Au. These 
values are in agreement with knowledge from solid state 
physics p3 24 as well as with the integration of optical 
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aed in detail in |16|]. As already known, the 
results obtained from the plasma model departs from the 
more accurate integration of optical data for small dis- 
tances. In this limit however, the thermal corrections do 
not play a significant role. In the present paper we will 
restrict our attention to the plasma model and discuss 
the validity of the results obtained in this manner at the 
end of the next section. 

We will also focus the attention on mirrors with a 
large optical thickness for which the reflection coefficients 
r±_ {iu, in) and rn {iu, in) correspond to a simple vacuum- 
metal interface. With the plasma model, these coeffi- 
cients are read as 
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For wavevectors ck smaller than ujp , mirrors may be con- 
sidered to be perfectly reflecting. When converted to the 
distance domain, this entails that the force approaches 
the ideal Casimir expression when evaluated at large dis- 
tances L 3> Ap. 

The Casimir energy will be obtained from the force by 
integration over the mirrors relative distance 



E 



F{x)dx 



(11) 



As this procedure is performed at constant temperature, 
the energy thus obtained corresponds to the thermody- 
namical definition of a free energy. For simplicity we will 
often use the denomination of an energy. We will define 
a factor rjp measuring the whole correction of energy due 
to conductivity and temperature effects with respect to 
the ideal Casimir energy 



E Ca 



E 

Ecas 

HcAtt 2 
720L 3 



(12) 



The positive value of the energy here means that the 
Casimir energy is a binding energy while the positive 
value of the force is associated with an attractive charac- 
ter. We will then define 2 factors ij^ and rj^ associated 
with each effect taken separately from each other, as in 

(4) . As already done for the force correction factors in 

( 5 ) , we will finally evaluate the quantity <5e which charac- 
terizes the degree of validity of the approximation where 
both effects are evaluated independently from each other. 
As mentioned in the Introduction, the results obtained 
for energy allows one to deal with the Casimir force in 
the plane-sphere geometry when trusting the proximity 
force theorem. 



III. NUMERICAL EVALUATIONS 

In the following we present the numerical evaluation 
of the correction factors of the Casimir force and energy 
using equations written in the former section. 

The force correction factor was evaluated for the exper- 
imentally relevant distance range of 0.1-10/^m with the 
help of equation (|7|) , supposing explicitly a plasma model 
for the dielectric function, and the result was normalized 
by the ideal Casimir force. A double integration over fre- 
quencies and wavevectors had to be performed. Due to 
the cosine dependence in (^|), the integrand turned out 



to be a highly oscillating function. Hence, the integra- 
tion required care although it was performed with stan- 
dard numerical routines. The energy correction factor 
was then calculated by numerically integrating the force 
and normalizing by the ideal Casimir energy (see equa- 
tion ([l2])). Integration was restricted to a finite interval, 
the upper limit exceeding at least by a factor of 10 4 the 
distance at which the energy value was calculated. Ex- 
tending the integration range by a factor of 100 changed 
the numerical result by less than 10~ 7 . 

The results of the numerical evaluation of rj-p are shown 
as the solid lines in figures [l] for Al and for Cu-Au as- 
suming a temperature ofT — 300AT. They are compared 
with the force reduction factor r]p due to finite conduc- 
tivity (dashed lines) and the force enhancement factor 
r)p calculated for perfect mirrors at 300K (dashed-dotted 
lines). 
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FIG. 1. Force correction factor for Al (upper figure) and 
Cu and Au (lower graph) as function of the mirrors distance 
at T = 300K. 

Figure || shows similar results for the factor rfE ob- 
tained through numerical evaluation of the Casimir free 
energy. The shape of the graphs is similar to the ones of 
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the force. However, while finite conductivity corrections 
are more important for the force, thermal effects have a 
larger influence on energy. 
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FIG. 2. Energy correction factor for Al (upper figure) and 
Cu and Au (lower graph) as function of the mirrors distance 
at T = 300^. 

For the force as well as for the energy, temperature cor- 
rections are negligible in the short distance limit while 
conductivity corrections may be ignored at large dis- 
tances. The whole correction factor rj behaves roughly 
as the product rfrj 1 of the 2 correction factors evaluated 
separately. However, both correction factors are appre- 
ciable in the distance range I — 4/im in between the two 
limiting cases. Since this range is important for the com- 
parison between experiments and theory, it is necessary 
to discuss in a more precise manner how good is the of- 
ten used approximation which identifies r\ to the product 
rfrj 1 -. In order to assess the quality of this approxima- 
tion, we have plotted in figure |^ the quantities 8-p and 
<$e as a function of the distance for Al, Cu-Au and two 
additional plasma wavelengths. A value of 8 — would 
signify that the approximation gives an exact estimation 
of the whole correction. An important outcome of our 



calculation is that the errors <5f and <5e are of the order 
of 1% for Al and Cu-Au at a temperature of 300iT. 
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FIG. 3. 8f (upper graph) and 6e (lower graph) as a func- 
tion of the mirrors distance. The results are given for the 
three metals Al, Cu-Au and two larger plasma wavelengths. 

For estimations at the 5% level, the separate calcula- 
tion of r] P and ry T and the evaluation of r\ as the product 
rfrj 1 can therefore be used. However, if a 1% level or 
a better accuracy is aimed at, this approximation is not 
sufficient. It should be noticed furthermore that the er- 
ror increases when the temperature or the plasma wave- 
length are increased. It becomes of the order of 4% for a 
plasma wavelength of 0.5 (im at 300K. The sign obtained 
for 8 means that the approximation gives too small values 
of force and energy. 

We want now to emphasize a few points. In order to 
make the discussion precise, we give numerical values of 
the correction factors for 2 experimentally relevant dis- 
tances, namely 0.5/im and 3/xm. The first distance cor- 
responds to the smallest distance for which the plasma 
model gives results in correct agreement with the integra- 
tion of optical data fll6|| . For this distance, the thermal 
corrections do not play a significant role (?yp = 1.000; 
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At shorter distances the results obtained with the plasma 
model depart from the values calculated from the integra- 
tion of optical data by more than 1%. Hence, the values 
of r]p and ?7g used for distances smaller than 0.5/zm have 
to take into account the more accurate dielectric function 
obtained through an integration of optical data |l(J . In 
this short distance range however, the whole correction 
factors 7/f and tje may be obtained as the products ?/p7yp 
and rferg. 

In the long distance range in contrast, the tempera- 
ture correction becomes predominant. The conductivity 
correction has still to be accounted for but it may be cal- 
culated by using the plasma model. This is illustrated 
by the correction factors obtained for a distance of 3/Ltm 
1.470). 
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For this distance, all corrections have to be taken into 
account. The metals cannot be considered as perfect re- 
flectors yet, the temperature corrections are significant 
and the deviation between the exact correction and the 
mere product has to be included if a high accuracy is 
aimed at. This is especially true in the case of Casimir 
energy. 



IV. SCALING LAWS FOR THE DEVIATIONS 

An inspection of figure |^ shows that the curves corre- 
sponding to different plasma wavelengths Ap have similar 
shapes with a maximum which is practically attained for 
the same distance between the mirrors. The amplitude 
of the deviations, which is larger for the energy than for 
the force, is found to vary linearly as a function of the 
plasma wavelength Ap . 

This scaling property is confirmed by figure ^ where we 
have drawn the deviations after an appropriate rescaling 



Ap 



The curves obtained for Ap and Ae for different plasma 
wavelengths at temperature T = 300-fT are nearly per- 
fectly identical to each other. These curves correspond 
to values of the plasma wavelength small compared to 
the thermal wavelength and the scaling law would not be 
obeyed so well otherwise. 
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FIG. 4. The deviations are represented for the force (up- 
per graph) and the free energy (lower graph) after the rescal- 
ing described by equation (15j). Different plasma wavelengths 
lead to nearly identical functions, drawn as dotted, dashed 
and dotted-dashed lines. These functions are hardly distin- 
guishable from the solid lines which represent the analytical 
expressions derived in the next section. 

In other words, the deviations Sf and §e are propor- 
tional to the factor ^ on one hand, and to the functions 
Ap and Ap on the other hand. The latter functions, 
which no longer depend on Ap , provide a simple method 
for reaching a good accuracy in the theoretical estimation 
of the whole correction factor 
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This method is less direct than the complete numerical 
integration of the forces which has been performed for ob- 



5 



taining the curves presented in the previous section. But 
it requires easier computations while nevertheless giving 
accurate estimations of the correction factors. Typically, 
the deviation S with a magnitude of the order of the % 
may be estimated with a much better precision through 
the mere inspection of figure ^. Alternatively, one may 
use the analytical expression of the functions A presented 
in the next section and drawn as the solid lines on figure 



V. ANALYTICAL EXPRESSIONS OF THE 
DEVIATIONS 

The results of numerical integrations presented in the 
foregoing section have shown that the deviations 5f and 
<5e are proportional to the plasma wavelength Ap, for 
plasma wavelengths small compared to the thermal wave- 
length. In this final section, we explain this scaling law 
by using a partial analytical integration of the whole cor- 
rection factors. 

To this aim, we write the force correction factor by 
dividing <^) by the value of the ideal Casimir force 



Vf = Vf + (Vf - !) + Ar/F 
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The first term in ( |17|) corresponds to the vacuum contri- 
bution (@) 
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with / still given by ([|). A dimensionless frequency 
y = — measured with respect to the wavevector has 
been introduced. Note also that the wavevector k is in- 
volved through the dimensionless quantity kL, except in 
the expressions of reflection coefficients. In (jig), the in- 
tegration over y may be performed analytically (see the 
appendix |X|). At long distances, ryp tends to the limit of 
perfect reflection with a known correction fl5f] 
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This expansion has been the subject of a number of pa- 
pers and it has been used to propose interpolation formu- 
las |2^,|2l[. However such a series expansion can hardly 
reproduce the behavior at small distances where rjp varies 
as which just means that the conductivity effect is 
not a small perturbation at short distances (see the ap- 
pendix |a|) . 

Coming back to the whole expression ( |l7j ) of the force 
correction factor, it remains to discuss the thermal con- 
tributions, that is the second and third terms. These two 
terms come from the contributions m ^ to (Q). The op- 
posite values of m give equal contributions and they have 
been gathered. The thermal contributions have been split 
in two parts, the second and third terms in (17), which 



correspond respectively to the limit of perfect mirrors on 
one hand 



240L 4 f 3 f 1 

~ / < \ ^ I dy cos (myKAT) f\ 

n m =i J o Jo 



■> 1 6 2kL _ I 

and the remainder on the other hand 



(20) 



240-L 4 °° f°° f 1 
ArjF = 1 — / dn k 3 dy cos (mynX^) A/ 



A/ = / - h 

„2kL 



+ 



e Ml. 



(21) 



The contribution (po| ) has been denoted (rjj — l) with 
?7p the correction factor obtained for perfect mirrors at a 
non zero temperature. For this term the integration over 
y is trivial and the integration over k may be performed 
analytically, leading to the known expression [ L3|— p5| 
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To obtain the overall correction factor ( |l7| ) it now re- 
mains to evaluate the last expression @. This can be 
done numerically, thus leading to the same results as in 
the previous section since no approximation has been per- 
formed up to now. But the results of the previous section 
suggest that we may obtain an accurate estimation of this 
term through an expansion in powers of Ap. The plasma 
wavelength Ap is indeed much smaller than the thermal 
wavelength At in all experimental situations studied up 
to now. Also, the deviation studied in the foregoing sec- 
tion is appreciable only for distances L much larger than 
Ap . Hence an accurate description of the deviation factor 
should be obtained by evaluating A77F at the first order 
in Ap. 

This first order term is easily deduced from (|lO|,pl"|) 
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It is proportional to Ap and to a function which does 
no longer depend on Ap 
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Collecting the results obtained up to now, we get an esti- 
mation of the force correction factor r]F valid in the long 
distance range L 3> Ap 
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Coming back to the notations of the previous section, 
this result is equivalent to the following expression for 
the function Ap 
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This function is plotted as the solid line on figure [| and it 
is found to fit well the results of the complete numerical 
integration presented in the previous section. 

Similar manipulations can be done for evaluating cor- 
rection factors for the Casimir free energy. We give be- 
low the main results, that is the thermal correction factor 
evaluated for perfect mirrors 
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and the first order correction 
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Since the long distance expansion of r/g up to first order 
in the plasma wavelength is given by 
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we deduce the function Ae 
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This function is plotted as the solid line on the second 
graph of figure |^ and also found to fit well the results of 
the numerical integration. 



VI. SUMMARY 

In the present paper, we have given an accurate evalua- 
tion of the Casimir force and Casimir free energy between 
2 plane metallic mirrors, taking into account conductiv- 
ity and temperature corrections at the same time. The 
whole corrections with respect to the ideal Casimir for- 
mulas, corresponding to perfect mirrors in vacuum, have 
been characterized by factors 7?p for the force and t?e for 
the energy. These factors have been computed through a 
numerical evaluation of the integral formulas. They have 
also been given a simplified form as a product of 3 terms, 
namely the reduction factor associated with conductivity 
at null temperature, the increase factor associated with 
temperature for perfect mirrors, and a further deviation 
factor measuring a kind of interplay between the two ef- 
fects. This last factor turns out to lie in the 1% range 
for metals used in the recent experiments performed at 
ambient temperature. Hence the conductivity and tem- 
perature corrections may be treated independently from 
each other and simply multiplied for theoretical estima- 
tions above this accuracy level. 

However, when accurate comparisons between exper- 
imental and theoretical values of the Casimir force are 
aimed at, the deviation factor has to be taken into ac- 
count in theoretical estimations. The deviation factor is 
appreciable for distances greater than the plasma wave- 
length Ap but smaller or of the order of the thermal wave- 
length At- We have used this property to derive a scaling 
law of the deviation factor. This law allows one to obtain 
a simple but accurate estimation of the Casimir force and 
free energy through a mere inspection of figure ^. Alter- 
natively one can use analytical expressions which have 
been obtained through a first order expansion in Ap of 
the thermal contributions to Casimir forces and fit well 
the results of complete numerical integration. 

We have represented the optical properties of metals 
by the plasma model. This model does not lead to re- 
liable estimations of the forces at small distances but 
this deficiency may be corrected by using the real di- 
electric function of the metals. This does not affect the 
discussion of the present paper, except for the fact that 
the pure conductivity effect has to be computed through 
an integration of optical data for distances smaller than 
0.5/im. Finally surface roughness corrections, which have 
not been considered in the present paper, are expected to 
play a significant role in theory-experiment comparisons 
in the short distance range. 
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APPENDIX A: THE VACUUM CONTRIBUTION 

In the present appendix, we give further analytical ex- 
pressions for the correction factor r/p due to conductivity, 
calculated with the plasma model for a null temperature. 

Introducing the notations 
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we rewrite the reflection coefficients (|To| ) 
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r± = -p 



r \\ = P 



In this case one integration may be performed analyti- 
cally in df) 



120L 4 f°° 3 2p 2 +pe KL g 

C1KK 
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+ a± 1 
arctan — 



+ 1 

— arctan — 

a i a i 



a± 



'e KL ± p 1 

b kL if pi 



At the large distance limit, rjp tends to unity, that is 
the value obtained for perfect reflectors. At the small 
distance limit, rjp is found to vary as Q 



L < A F 



[10 

[11 

[12 

[13 
[14 
[15 

[ie; 
[17; 
[is; 

[19 

[20 

[21 

[22' 
[23 
[24 



E. Fischbach and C. Talmadge, The Search for Non New- 
tonian Gravity (AIP Press/Springer Verlag, 1998) and 
references therein; E. Fischbach and D.E. Krause, Phys. 
Rev. Lett. 82 4753 (1999) 

G. Carugno, Z. Fontana, R. Onofrio and C. Rizzo, Phys. 
Rev. D55 6591 (1997) 

M. Bordag, B. Geyer, G.L. Klimchitskaya, and V.M. 

Mostepanenko, Phys. Rev. D60 055004 (1999) 

E.M. Lifshitz, Sov. Phys. JETP 2 73 (1956); E.M. Lif- 

shitz and L.P. Pitaevskii, Landau and Lifshitz Course of 

Theoretical Physics: Statistical Physics Part 2 ch VIII 

(Butterworth-Heinemann, 1980) 

J. Mehra, Physica 57 147 (1967) 

L.S. Brown and G.J. Maclay, Phys. Rev. 184 1272 (1969) 
J. Schwinger, L.L. de Raad Jr., and K.A. Milton, Annals 
of Physics 115 1 (1978) 

A. Lam brecht and S. Re ynaud, Eur. Phys. J. D, to ap- 



pear; yiant-ph/9907105 



Lifshitz, Quart. 

Tsang, 



B.V. Deriagin, LI. Abrikosova and E.M. 
Rev. 10, 295 (1968) 

J. Blocki, J. Randrup, W.J. Swiatecki and C.F. 
Ann. Physics 105, 427 (1977) 

V.M. Mostepanenko and N.N. Trunov Sov. J. Nucl. Phys. 
42 812 (1985) 

V.B. Bezerra, G.L. Klimchitskaya and C. Romero, Mod. 
Phys. Lett. 12, 2613 (1997) 

G. L. Klimchitskaya, A. Roy, U. Mohideen and V.M. 
Mostepanenko, Phys. Rev. A60 3487 (1999) 

M.T. Jaekel and S. Reynaud, J. Physique 1-1 1395 (1991) 
L.G. Schulz, Phil. Mag. Suppl. 6 102 (1957) 

H. Ehrenreich and H.R. Philipp, Phys. Rev. 128 1622 
(1962) 



7? F ~ a— 
Ap 



30 



1.193 



dK e 



/sinh f 



A' 2 



cosh y 



[1] H.B.G. Casimir, Proc. Kon. Nederl. Akad. Wet. 51 793 
(1948) 

[2] B.V. Deriagin and I.I. Abrikosova, Soviet Physics JETP 
3 819 (1957) 

[3] M.J. Sparnaay, Physica XXIV 751 (1958); W. Black, 
J.G.V. De Jongh, J.Th.G. Overbeek and M.J. Sparnaay, 
Transactions of the Faraday Society 56 1597 (1960) 
[4] D. Tabor and R.H.S. Winterton, Nature 219 1120 (1968) 
[5] E.S. Sabisky and CH. Anderson, Phys. Rev. A7 790 
(1973) 

[6] S.K. Lamoreaux, Phys. Rev. Lett. 78 5 (1997); erratum 

in Phys. Rev. Lett. 81 5475 (1998) 
[7] U. Mohideen and A. Roy, Phys. Rev. Lett. 81 4549 (1998) 
[8] A. Roy, C. Lin and U. Mohideen, Phys. Rev. D 60, 

111101 (1999) 



8 



